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ABSTRACT 


^->The  research  4es«r1vrri  in  -hie  repeat  ieals  vrith  the  application  of 
abstract  techniques  in  the  3tudy  of  the  behavior  and  asymptotic)  stability 
of  solutions  to  differential  equations  of  reaction-diffusion  type.  These 
equations  include  models  in  a  cellular  control  process  with  positive  feed¬ 
back  and  a  model  in  *as-liquid  reactions.  Most  of  the  results -aim  based  on 
techniques  that  involve  the  determination  of  invariant  sets  for  Solutions  and 
the  preservation  of  inequalities  between  solutions. 


1  JAN  73 


1473  EDITION  OF  1  ITOV  65 
IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


5 


where  a,  >  0  Tor  i  *  1 . ,m,  <J>1,  and 

(2)  h  (r)  *  r^d+r0)”1  for  all  r>0. 

0  “* 

Several  authors  have  investigated  the  existence  of  and  the  stability  properties 
of  ncnnegative  critical  points  (equilibrium  solutions)  to  this  feedback  system. 
However,  the  earlier  works  on  this  model  did  not  use  the  fact  that  solutions  to 
this  equation  preserve  inequalities:  if  z(.;n)  is  the  solution  to  (1)  for  each 

•jj 

given  nonnegative  initial  value  n  «  (n,)^  »  then  z(t;n)  >_  z(t;TT)  (componentwise) 

for  all  t  >  0  whenver  (componentwise).  This  observation  allows  one  to 

both  improve  and  simplify  previous  results.  Also,  abstract  techniques  allow  much 
more  general  equations  to  be  studied.  First  we  consider  the  associated  reaction- 
diffusion  system. 

’•A  '  VxA  -  Vi  *  h«(“a) 

3.^  -  dgJ^Ug  “  a2U2  *  'h  (***)  C  tO,-)x(0,l) 

(3)  . 

3  u  *  i  3  u  -a  u  ♦  u  , 

t  m  m  xx  m  mm  m-1 

subject  to  the  initial  conditions 

(1»)  u^Q.x)  ■  x,  (x)  >_0, . . .  ,ua(0,x)  ■  X, 'x)>p  for  xc(0,l) 
and  the  boundary  conditions 

po3xUi^t’°*  -  qQU1(t,0)  *  p13jtu1(t,l)  *  q^U.D  -  0 

(5) 

for  all  t>0  and  i  *  l,...,m 

where  qQ,  q,  >_0  and  p.  ,p^  c  (0,  1}  are  such  that  q^  *  1  whenever  p  ■  0  and 
q,  *  1  whenever  p,  ■  0.  Consider  also  the  linear  eigenvalue  problem 

ft 

*  (x)  ♦  \t(x)  ■  0,  0<x<l 

(6) 

P0M°)  -  q^HO)  •  pV(l)  ♦  qx»(l)  *  0 


and  let  X  *  X^  denote  the  first  eigenvalue  (note  that  X^  >_  0  and  that  X^  *  0 

'  i 

only  if  the  boundary  conditions  are  ♦  (0)  *  ♦  .1)  -  0).  The  techniques  and 
results  in  (3)  show  that  the  following  behavior  pattens  for  solutions  to  (3)  - 
(5)  in  valid: 

THEORIW  1.  Suppose  that  o *  1  in  the  definition  (2)  of  h  .  that  d_  . d  are 

““ ~ ~  o  i  a 

positive  constants  in  equation  (3),  and  that  X^  i3  the  first  eigenvalue  of 
equation  (6). 


(i)  If  n  (a.  ♦  X  d,  )  1  then  for  each  nonnegative  initial  value 

i-1  1  1  1 

X  *  the  solution  u  -  (u^)25,  to  (3)  -  (5)  exists  on 

(0,«)  x  (0,11  and  u(t,x)  ■*  0  a3  t  ■»  »  ’uniformly  for  xe  [0,1]. 

(ii)  If  “  (a.  ♦  X,d()  <  1  then  there  is  a  ’unique  nontrivial,  nonnegative 

i-1 

equilibrium  solution  *  -  (i,)31,  to  (3)  -  (5): 

0  -  d-*,  (x)  -  a.*,(x)  ♦  h.(*  ( x ) ) 

xi  Xl  X  a 

0  *  &2*2  ^  ~  ®2*2^x^  *  0<x<l 


d  * 
a  a 


(xl 


0(0 

21  12 


*  .  (x) 

a-1 


and 

9o*\(0)  -  <1^(0)  *  pi* '^(1)  ♦  q.*,(ll  *  0  i-1, . . .  ,a 

Moreover,  *^(x)  >  0  for  ail  x  and  i,  and  fbr  each  ncnnegative , 
nontrivial  initial  ’value  x  ■  (x^-0,  the  solution  u  -  (u4)m.  exists 
on  (0,-)  x  (0,1)  and  u(t,x)  *  (x)  as  t  •*  •,  uniformly  for  xe(0,l]. 


It  is  of  interest  to  note  that  when  o  -  1  the  behavior  of  the  solutions  to  the 
ordinary  differential  equation  (1)  is  precisely  the  same  as  that  of  the  reaction- 


r 


i 
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diffusion  model  if  one  replaces  a,  by  a,  ♦  A^d^  in  (1). 

A  result  similar  to  Theorem  1  can  be  obtained  by  considering  time  delays  in 
the  feedback  terms.  Again  suppose  o*  1  and  consider  the  delay  differential 
equation. 


(7) 


'(t) 


-  -  Vl(t)  *  hl(VmCt-rm)K 

■  -  «■  v,  (t-r^) , 


w  (.)  *  ^(a),  - r^ <s<0 

w  (.)  *  <^(s),  - r^<s<0 


•  *  (t) 

3  3 


*  Vl(t"r«-1)' 


w  (s) 

3 


■r  <s<0 

3 - 


where  r.  >_  0  and  i.  is  a  giver,  nonnegative  function  on  [-r.  ,0]  for  each  i  *  l,...,a 

TV.e  results  for  the  system  (7)  are  completely  analogous  to  the  corresponding 
reaction-diffusion  system  (3)  -  (5): 

THE: REM  2:  Consider  the  time  delay  differential  equation  (7). 

i)  If  1'  i,  >  1  then  for  each  nonnegative  initial  function  t  *  (  e,  )B^ 

i-1 

the  solution  w  ■  v. )a,  to  (7)  exist  on  [0,»)  and  w(t)  *  9  as  t 

(ii  If  B  a  <  1  then  there  is  a  unique  nonnegative,  nontrivial  vector 
*»*  1 

i-1 

6  -  { 6.  )a  such  that  w(t)  35  for  t>0  (and  *.  (s)  36.  on  [-r,  ,0] 

"“111 

for  i  -  1, ...,»)  is  a  constant  solution  to  (7).  Moreover,  6^  >  0 

for  all  i  -  l,...,m  and  if  w  is  a  nontrivial  solution  to  (7)  then 
v(t)  -*■  5*  as  t  ■*  -.  This  result  is  from  the  article  (6). 

Results  have  also  been  obtained  on  the  behavior  of  solution  to  a  nonlinear 
parabolic  system  that  arises  in  the  study  of  gas-liquil  reactions.  The  fundamental 


equation  has  the  form 
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(8) 


*tu(t,x)  •  <^3  u(t,x)  ~  kY(u(t,x),  v(t,x)) 


t  >  0 


3tv(t,x)  ■  i.,3  v(t,x)  ♦  k  y  (u(t ,x) ,  v(t,x))  0 <x<o 


subject  to  the  boundary  conditions 


(9) 


u(t,0)  *  a 


3  v(t,0)  *  0 


O  X 

u(t,o)  -  0  v(t,o)  »  0 

and  the  Initial  conditions 


t  >  0 


(10)  u(0,x)  »  u.  ( x ) ,  v(0,x)  »  v  (x), 

O  0 


0  <x<o. 


Here  it  is  assumed  that  a^ ,  bQ  and  k  are  given  positive  constants,  c  is 


a  positive  number  or  ♦  •,  and  y  is  a  real  value  C*-  function  on  (0,  a^] 


either 
x  (0,  bQ) 


such  that 

'y1)  y(0,v)  *  Y(u,b  )  *  0  for  all  (u,v)  c  [0,  a  ]  x  (0,  b  ]  and 
o  oo 

(y2)  3  y(u,v)  >  0  and  3^y(u,v)  <  0  for  all  (u,v)  c  (0,  a^  x  (0,  bo). 

The  typical  example  of  such  a  function  Y  is  y(u,v)  ■  vi  ( b.-v)1  where  p,q  >_  l). 

The  initial  conditions  u  and  v  are  al30  assuaed  to  satisfy  0  <  u  (x)  <  a  and 

oo  —  o  —  o 

0  <  v  (x)  <  b  for  alacst  all  x  c  [0,o]. 

—  o  —  o 

The  analysis  of  this  3y3ten  fits  nicely  into  an  abstract  setting,  let  Lx  * 

1  2 

I/K 0, a)  i  K"  1  be  the  Banach  space  of  all  measurable  functions  f*  t,  ,  ,  from 

2 

(0,  a)  into  H  3uch  that 

I  Ul  Ix  5  dj  (!  *1(x)|-*‘  I  *2(x)|  )  lx  <  - 

and  define  3  C  L*  by  D  *  ((♦,,  ♦,)  t  L^:  0  <  *  (x)  <  a  and  0  <  ♦-(x)  <  b  a.e. 

It  “O  —  d  *  o 

on  (0,  o)}.  Note  that  D  is  a  closed  convex  subset  of  Lx  (with  empty  interior!) 
and  that  D  is  bounded  only  in  case  o  <  •. If  (♦,  ,  ft,)  and  (*,  ,*0)  are  in  we 
write  *  <_  only  in  case  *^(x)  <_  *  (x)  and  *^(x)  <_  a.e.  on  (0,o).  Concerning 

the  fundamental  existence  and  behavior  of  solutions  to  (8)  -  (10)  we  have  the 


following  result : 


THEOREM  3.  Suppose  that  (yl)  and  (y2)  are  satisfied.  Then  for  each  initial 

value  (u  ,  v  )  e  D  the  system  (9)  -  (10)  has  a  unique  solution  (u,v)  and  this 
o  o 

solution  exists  on  (0,  •  )  x  (0,<j)  and  satisfies  (u(t,.),  v(t,.))  e  D  for  all 

t>0  (i.e.,  0<u(t,x)<a  and  0<v(t,x)<b  for  all  t>0  and  0<x<c).  Moreover, 

—  — o  —  -  o 

i f  S(t)(u  ,v  )  5  ( u( t,.),v(t,.))  for  each  initial  value  (u  ,v  )  in  D,  then  S  * 
o  O  0  0 

{S(t):t>_0}  is  a  semigroup  of  nonlinear  operators  mapping  D  into  D  that  is  non- 

expansive  and  order  preserving: 

(i)  S(t-*-s)(u  ,v  )  »  S(t)S(s)(u  , v  )  for  all  t;3>0  and  (u  ,v  )  c  D; 

0  0  0  0  -  0  0 

;ii;  !  ]S(t)(uo,v3)  -  S(t)(u1,v1)  I  '  I  (uQ,v0)  -  (u^vJll,  for  all 

t>0  and  (u  ,v  ),  (o,,v .  )  e  D;  and 

( iii)  S(t)(u  ,v  )  <  S(t)(u,  ,v  )  whenever  t  >  0  and  (u  ,v  ) ,  (u.  ,v,  )  e  D  with 

00”  11  ”  00  11 

(u  ,V  )  <  (u.  ,V,  )  . 

0  0”  1  I 

The  proof  of  this  theorem  follows  from  the  results  in  [8]  and  i3  indicated  in  the 
article  [5). 

The  behavior  of  solutions  to  (8)  -  (10)  a3  t  -*•  •  has  also  be  investigated  in 
(5)  as  well  as  in  (1).  The  analysis  of  this  behavior  divides  naturally  into  two 

parts:  The  film  theory  (when  a  <  •)  and  the  penetration  theory  (when  a  *  •). 

THEOREM  1* .  Suppose  that  the  conditions  and  notations  in  Theorem  3  are  satisfied. 

(i)  If  a  <  m  there  is  a  unique  equilibrium  solution  *  ■  (♦. .  *„)  to  (8)  - 

(10)  (i.e.,  1,  <>.  -  ky(*)  *  0  and  *  *  Y  (♦)  *  0  on  [0,o], 

(0)  *  a  ,  (0)  *  0,  and  (o)  *  *~(a)  *  0]  and  for  each  (u  ,v  )cD, 

10  c  1  d  0  o 

S(t)(u3,v^)  ♦  ms  t  •  uniformly  for  xc  (0,o). 

(ii)  If  a  ■  •  there  is  no  equilibrium  solution  to  (8)  -  (10)  and  for  each 

(u  ,v  )  t  D,  (S(t)(u  ,v  ))  (x)  ■*  (  a  ,b  )  as  t  +  •  uniformly  for  x  in  each 

00  00  0  0 


bounded  subset  of  (0,  •). 


The  results  in  Theorem  4  as  well  as  more  general  results  can  be  found  in  the 
paper  (5). 
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In  (2]  abstract  techniques  are  applied  to  a  mathematical  model  of  a  gas 
exchange  system.  This  system  can  be  put  in  the  following  form.  Suppose  that 
o ,  t,  P,  o,,  and  k,  ,  ,  k^,  are  positive  numbers;  that  S  ■  (S^3.,  and 

y  ■  ( y  )3  are  funct ions  from  3?'  into  3F';  and  that  u  ■  (u^)3  and 

w  *  (v^)3.  are  C3,  on[D,«) .  If  zaz,  )3^  and  w  »  (w,  ",  map  [0,1]  into  1R',  the 

problem  i3  to  determine  the  existence  of  functions  u  *  (u^)3,  and  v  =  (v^)3^ 

from  (0,  •  ] x  [0,  1]  into  JR3  and  a  function  c  from  (0,«*)  x[0,l]  into  I  such  that 


t  >  0,  0<x<l 
and  i  *  1,  2,  3. 


is  satisfied  along  with  the  boundary  conditions 

u(t,0)  ■  u(t),  v(t ,0)  *  w ( t )  n 

(12) 

3xu(t,l)  *  3^v( t ,x)  *  c(t,l)  *  0 
the  initial  conditions 

(13)  u(0,x)  ■  x(x)  and  v(0,x)  *  v(x)  0<x<t, 

as  well  as  the  side  condition 

(1U)  u^( t; *x)  ♦  Ujft.x)  ♦  ujt.x)  5  P  for  t>0,  0<x<t. 

The  following  existence  result  can  be  found  in  [2]: 


In  addition  to  the  suppostions  and  notations  in  the  preceding  paragraph, 


let  S,  ,  P„,  and  be  positive  numbers,  let 

Ax  -  u  c  S3:  ♦  Cj  ♦  C3  -  P  and  ^  *3  >0) 

A2  -  U  c  E3:  0«  C1iR1  for  i  »  1,  2,  3), 


li 


and  assume  also  that,  the  following  is  valid: 

(HI)  (z(x),  v(x))e  A, xA,  for  0<x<l  and  (u(t),  v(t))e  A,  xA,  for  all 

t  >.0;  and 

(H2)  If  U.n)  e  AxxA^  and  J  e  (1,  2,  3)  then  5,  *  0  implies  BjU.n) 

n  »  0  implies  Y,(£,n)>.0;  and  n,  *  R,  implies  y,(C,  n/f.0. 

Then  the  system  (10)  -  (lM  has  solution  on  [0,«)x  [0,1]. 

As  opposed  to  considering  this  system  as  a  differential  equation  with  side 
condition,  it  is  shown  in  [2]  that  one  can  establish  the  existence  of  solutions  by 
•using  the  abstract  theory  of  invariant  sets  for  evolution  systems  (see  [9]). 

In  [U]  the  [strict,  componentwise)  positiveness  of  reaction  diffusion  systems 
of  the  form 

♦ 

■ .  j.  .  J.u,  *  d,  du,  *  ?.  u  ,  Vu,  '  on  (0,*;x  Q  for  i  *  1 , . . .  ,m 

with  initial  and  boundary  conditions 

(16)  u  *  9  on  (0,»)x3fi  and  u  *  x  on  (0)  x  fl 

are  considered  as  well  as  extension  of  such  systems)  where  fl  is  a  smooth  bounded 
iomain  in  3  ,  A  is  the  Laplacian  operator  on  Q,  7  is  the  gradient  operator,  and 
d,  is  a  positive  constant  for  each  i  *  l,...,m.  Criteria  is  given  to  insure  that 
if  the  initial  value  x  is  nontrivial  componentwise  nonnegative  on  then  each 
component  of  the  solution  u  is  strictly  positive  for  positive  time:  u  (o,x)>_  0 
but  not  identically  zero)  implies  that  u^!t,x)»0  for  all  t>0,  xefl  and  i  *  l,...m. 

Thi3  results  are  obtain  by  comparison  with  the  ordinary  differential  equation. 

(17)  z  ’(t)  ■  g(t,z(t)),  t  >  0 

where  g  *  (g^)a,  is  defined  on  [0,«)x  J1  by 

g.(t,C)  *  F  (t,x  9)  for  all  (t,£)  c  [O/lxS31  and  i  ■  l,...,m. 

X  1  o 

(the  point  xq  is  some  given  point  in  fi).  ITiis  type  of  maximum  principle  extends 
a  similar  result  established  in  [ 3 1 • 
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